Dynamical behavior of disordered spring networks 
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We study the dynamical rheology of spring networks with a percolation model constructed by 
bond dilution in a two-dimensional triangular lattice. Hydrodynamic interactions are implemented 
by a Stokesian viscous coupling between the network nodes and a uniformly deforming liquid. 
Our simulations show that in a critical connectivity regime, these systems display weak power law 
rheology in which the complex shear modulus scales with frequency as G* ~ (i^)^ where A = 0.41, 
in discord with a mean field prediction of A = 1/2. The weak power law rheology in the critical 
regime can be understood from a simple scaling relation between the macroscopic rheology and the 
nonaffine strain fluctuations, which diverge with vanishing frequency for isostatic networks. We 
construct a dynamic effective medium theory, which describes the rheology at the mean field level 
and correctly predicts the scaling form for the rheology of networks with near-isostatic connectivities. 
Surprisingly, by connecting this critical scaling of the rheology with that of the strain fluctuations, 
we find that the dynamical behavior of disordered spring networks is fully determined by the critical 
exponents that govern the behavior of elastic network in the absence of viscous interactions. 



Disordered mechanical networks are used to model 
a variety of systems including network glasses [IHZ]j 
jammed packings [SHllj and semiflexible biopolymer net- 
works [HHTT] . Though the specific responses of different 
classes of networks may vary, the general character of 
the mechanical behavior depends on the networks' con- 
nectivity. Above the so-called isostatic connectivity, net- 
works are mechanically stable and resist static shear de- 
formations [IB]. The lack of mechanical rigidty below 
this connectivity is due to zero-energy, floppy deforma- 
tion modes; however, the mechanical response can be sta- 
bilized by additional weak interactions [3, 4, [TOl HSl [19] 
or internal stresses [201 [21] ■ Recently, focus has shifted 
towards the dynamic behavior of networks at the verge of 
mechanical stability with additional viscous interactions 
[221 - [25] . Marginally stable spring networks in jammed 
configurations exhibit a rich dynamic mechanical re- 
sponse in the presence of damping forces [25H?7] : in the 
vicinity of the isostatic connectivity, the dynamic shear 
modulus was found numerically to scale with frequency as 
•\/aJ, in accord with a mean field prediction. In contrast, 
bond-diluted lattice-based networks — a major class of 
mechanical systems distinguishable from jammed pack- 
ings [iT] — do not exhibit mean-field behavior [31-151 116j. 
This class of models has provided insight in both the lin- 
ear and nonlinear elastic response of disordered spring 
and fiber networks. However, little is understood about 
the dynamical response of such systems in the presence 
of viscous coupling between the network and a liquid. 

Here we study the linear rheological behavior of perco- 
lating disordered spring networks immersed in a Newto- 
nian liquid. Disordered spring networks are constructed 
on a triangular lattice in 2D, diluting bonds to control the 
networks' connectivity, ranging from well below to well 
above the isostatic connectivity. Analysis of the simu- 




FIG. 1. Exemplary samples of networks undergoing dynamic 
deformation, where the coloring of the plots is a qualitative 
representation of the nonaffinity of the two nodes connected 
by each bond (redder is more nonaffine) . Networks with p ^ 
Pc (ad), p ~ pc (be), and p <^ Pc (cf) are illustrated; the 
top row depicts high-frequency oscillation, uj ^ u* , and the 
bottom row depicts low-frequency oscillation, oj ^ lj*. In the 
high frequency limit, all networks deform affinely because the 
spring network couples strongly with the affine deformation 
of the viscous network (abc). In the low frequency limit, 
the spring network deformation is nonaffine, and the degree 
of nonaffinity increases as p approaches the isostatic bond 
probabilty Pc (def). 



lated complex shear modulus of such networks reveals 
three dynamical regimes. For connectivities well above 
isostaticity, the networks behave as a solid with weak 
viscous coupling to the fluid. The deformation of these 
networks is nonaffine at low frequencies, but becomes 
increasingly affine at high frequencies, as illustrated in 
Fig. [T] The second regime is in the vicinity of the iso- 
static connectivity, where we observe weak power law rhe- 
ology over a frequency range extending to zero- frequency 
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at the isostatic connectivity; this critical slowing down 
indicates diverging relaxation timescales. In this criti- 
cal regime the complex shear modulus scales with fre- 
quency as G* ^ («^)^, where A 0.41. In the third 
regime, with connectivities well below isostaticity, the 
networks behave as Maxwell fluids, crossing over from 
a fluid-like to a solid-like response at a connectivity- 
independent characteristic frequency. Our results are 
qualitatively consistent with results by Tighe p5H?7] on 
disordered spring networks based on jammed configura- 
tions, although in those systems the dynamic exponent 
A = 1/2 was found at near-isostatic connectivities. In 
addition, we develop a dynamic effective medium the- 
ory (EMT) for the rheological response of bond-diluted 
lattice-based networks which quantitatiely describes the 
regimes on both sides of the isostatic connectivity. This 
dynamic EMT serves as a framework that can be ex- 
panded to bond-bending and fiber networks O [6l ITSl - 
[iTl [19] . Our EMT calculation, taken together with scal- 
ing arguments, indicates that the dynamical properties 
of the networks are directly implied from the scaling of 
the strain fluctuations and mechanical behavior of purely 
elastic spring networks. 

The mechanical response of spring networks depends 
sensitively on the network's coordination number z, the 
average number of springs attached at a node in the 
network, not including dangling springs. Maxwell's con- 
straint counting argument indicates a critical condition 
of Zc — 2d for the onset of mechanical rigidity in a central 
force network in d dimensions [18] . To create a network 
with variable z in a range spanning from well below z^. 
to well above, springs are arranged on a triangular lat- 
tice and are removed with a probability 1 — p such that 
the network connectivity is roughly z ~ 6p, resulting in 
a network architecture distinguishable from jammed net- 
works [TT]. Using units in which the spring rest length 
^0 and stiffness /i are both 1, the energy can be written 
for small relative deformations Uij = Uj — between 
neighboring nodes i and j as 

^=\Y.a^M^r^^of (1) 

where gij = fi = 1 for a present bond or for an absent 
bond and f ^ is a unit vector directed along the ij-bond 
in the undeformed reference lattice. This network is em- 
bedded in a viscous fluid in the low-Reynolds number 
limit. As the network is deformed, hydrodynamic inter- 
actions between nodes are ignored, and the fluid deforms 
affinely; consequently, the net force on a node i is given 
by the viscous Stokes drag and the elastic central forces 
due to the springs to which it is connected, 

ii = 4:ni]a{ui - Ufiuid) + ^ giji^ij ■ (2) 

where lifluid is the velocity field of the underlying fluid, tj 
is the fluid's viscosity, and the summation of (j) is over 
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FIG. 2. Plot of a typical strain (blue) applied to the spring 
network and the resultant stress (red) of the network for 
p — 0.6 (a). We observe a phase shift S between the applied 
strain and the response, as well as a decrease in amplitude 
CTo/70, from which we can calculate the complex shear modu- 
lus given in Eq. Q. The effect of the dynamical shear strain 
is shown for a sample network on the right; (b) shows the ini- 
tial network state at to = 0, and (c) shows the network state 
after a quarter oscillation at t = tt/2uj. 

nearest neighbors of node i. For the 2D drag coefficient 
in the Stokes term, we associated a hydrodynamic radius 
a to a network node and chose natural units such that 
rja — 1. 

To study the rheology of these diluted networks nu- 
merically, we impose a time-dependent oscillatory shear 
strain along the two parallel sides of the network with fre- 
quency uj (Fig. [2]bc), which is applied by using periodic, 
Lees-Edwards boundary conditions [21]. The resulting 
macroscopic shear stress is calculated by 

where A is the surface area of the network and the 
x-component of the force between two nodes fij^x = 
gij{uij ■ rij)fij j.. A typical stress-strain relationship for 
an oscillating network is depicted in Fig. [2] a. From this, 
we determine the complex shear modulus 

G*{uj) = G'{Lj)+iG"{Lu) = — [cos{S{uj))+ism{S{oj))] , 

70 

(4) 

where ctq is the magnitude of the observed shear stress, 70 
that of the imposed shear strain, and 5{oj) the phase lag 
between the stress and the strain at frequency uj. Solid- 
like systems are dominated by G', the storage modulus, 
whereas liquid- like systems are dominated by G", the loss 
modulus. 

In the quasistatic limit a; — 0, the behavior of the net- 
work is determined by the relation between z and Zc or, 
equivalently, between p and Pc] the elastic shear mod- 
ulus vanishes continuously at the isostatic bond prob- 
ability Pc as G' ^ Ap-'', where the rigidity exponent 
/ ~ 1.4 [nifin] and Ap = p — pc- By contrast, in the high 
frequency limit w — cx), the network's response becomes 
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FIG. 3. Simulated rheology for various values of p and uj. All simulations were performed on a network of size 100 x 100. The 
frequency and storage (G') and loss (G") moduli are in natural units, as described in the main text. A slope of 2 in a) and 1 
in b) indicates low-frequency fluid-like behavior. At frequencies below u* near pc ~ 0.649, the shear modulus scales as a power 
law with exponent 0.41. Additionally, the EMT results for the storage modulus c) and for the loss modulus d) are shown, in 
reasonable agreement with the simulations. The insets show the low- and high-frequency limits according to Eqs. (171 and 
(l20l. 



affine, G — >■ Gaffino = P"\/3/4, for all bond probabilities. 
(Fig. [3] a). At high frequencies, nonaffine deformations 
are suppressed by large drag forces between the network 
and the affinely deforming fluid; thus, at high strain rates 
the fluid effectively dictates the behavior of the network. 

Networks with connectivities well below the isostatic 
connectivity behave as Maxwell fluids — the storage mod- 
ulus vanishes as uP' and the loss modulus as lo at low fre- 
quencies, crossing over at a frequency w* — l/47r, set 
by comparing the stretch modulus of a spring to the 
drag coefficient of a node, to high-frequency afline elastic 
behavior. By contrast, hyperstatic networks cross over 
at this characteristic frequency from a nonaffine to an 
affine solid-like gel at uj* . This transition is accompa- 
nied by a maximum in the loss modulus, G" (Fig.|3]b). 
At connectivities near isostaticity, the shear modulus 
appears to exhibit a power law regime G* ~ (iw)'^, 
where A « 0.41, extending to the zero-frequency limit as 



z — > Zc. Spring networks in jammed configurations dis- 
play a similar rheological behavior but with mean field 
exponent A = 1/2 [25H27] . The distinct regimes observed 
here over a connectivity range 3 < z < 6 are visualized 
by plotting the inverse loss tangent G'/G" for a range of 
network connectives and frequencies, as shown in Fig.[5]a. 

To supplement the complex shear modulus as a de- 
scription of the macroscopic behavior of spring networks, 
we study the fluctuations in the microscopic network de- 
formations using a simple one-point nonaffinity measure 
based on the nonaffine component of a node's displace- 
ment, 5\i: 




where the brackets indicate an average over network 
nodes and time. Studies on elastic networks have shown 
that the nonaffine fiuctuations diverge aX p = Pc like 



FIG. 4. The nonafRnity measure obtained from simulations as 
a function of frequency (main plot), collapsed by the scaling 

form r = |App^ *± (^cjjApp'^'^j (lower-left inset), and G" 

against Tuj (upper-right inset). At lj* — l/in, the nonaffinity 
crosses over from a high-frequency regime, in which nonafiine 
fluctuations are minimized by the fluid network with scaling 
exponent —2, to a low-frequency regime, where the nonaffinity 
measure scales with scaling exponent —S = —0.59. The extent 
of this latter region extends to zero frequency as the bond 
probability approaches the isostatic point. 



r = r±|Ap|-^ [TOl Hg Eg, where A w 2.2 for bond- 
diluted 2D triangular networks. In dynamical networks, 
nonaffine fluctuations will be suppressed by the affinely 
deforming viscous fluid. Consistent with prior work |26j . 
our simulations show that the nonaffine fluctuations ex- 
hibit a frequency dependence in certain regimes, as shown 
in Figs.|4]and[5]b: far away frompc, the nonaffine fluctu- 
ations are frequency-independent at low w. By contrast, 
near isostaticity, the nonaffinity measure depends on fre- 
quency as a power law F ^ w^'', with 6 = 0.59. For all 
values of p, the nonaffinity vanishes as F '-^ uj~'^ beyond 
the crossover frequency uj* . The frequency and connec- 
tivity dependence of the nonaffinity in dynamic networks 
can be captured by the scaling ansatz |26| : 

F = \Ap\-^ ^± (c./we) , (6) 

where the critical relaxation frequency, uic — lApj*^", de- 
scribes the slowest relaxation rate in the system. In- 
deed, we find a good collapse with this scaling form using 
A — 2.2 as determined previously in elastic networks |16| 
and (j)^ = 3.6, as shown in the lower inset of Fig. |4] 
Near isostaticity, the nonaffine fluctuations are finite, and 
F ^ w^*; therefore, it becomes clear that \I'±(a;) ^ x^^ 
with S — to eliminate the Ap dependence in Eq. 

The nonaffine fluctuations can be related to the shear 
modulus by estimating the dissipated power in the sys- 
tem in the critical regime in two different ways [261 12 7j . 
The viscous forces scale as Vt^o^^, giving rise to a dissi- 
pated power W ~ F(7oa;)^, while on a macroscopic level 
the dissipated power is given by = ^G"u}^q. It follows 
that G" ~ Tuj. To test this relation over a broad range 
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FIG. 5. The dynamical behavior of bond-diluted spring net- 
works. In the high frequency limit, lj > uj* , the network 
deforms affinely for all connectivities; therefore, G' 3> G" 
(panel a) , and the nonaffinity measure of the network F van- 
ishes as uj~^ (panel b). Below the crossover frequency uj* , 
networks with p ^ Pc are non-affine, solid-like gels with 
G' ^ G" , and networks with p <^ Pc exhibit a more fluid-like 
response with G' <C G" . In the neighborhood of pc, there is a 
critical regime where G' ~ G" ~ (iuj)^ where A = 0.41. The 
extent of this regime is bound by a|j|^| Ap]'*'^ , where a± has 
been inferred from the collapse described in Fig. |6] Panel b) 
shows log(r) for various p and uj, illustrating the divergence 
of the nonaffinity measure at pc in the low-frequency limit. 
The inset shows F on a linear scale. 



of connectivities and frequencies, we plot G" against Tuj 
and find that all data collapse onto a curve with a slope 
of 1, affirming this correlation (Fig. |4] upper inset). Fur- 
thermore, since near isostaticity G" ^ uj^ with A « 0.41 
(Fig. |3]ab), this connection between G" and F implies 
that the dynamical exponents 5 and A are related as 

<5 = 1 - A, (7) 

which is consistent with our observation oi 5 — 0.59 
(Fig. |4] inset) . 

DYNAMIC EFFECTIVE MEDIUM THEORY 

To provide insight into the simulated dynamic rhe- 
ology, we develop a dynamic Effective Medium Theory 
(EMT) m H Ha [ni nil EH [30]. THs approach is based 
on the construction of a mapping from a lattice network 
where the spring between nodes i and j has a spring con- 
stant gij, drawn from a probability distribution P(gij), 
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onto a lattice with uniform, frequency-dependent bond 
stiffness g(a;); this effective lattice mimics the mechani- 
cal response of the disordered network at the same global 
strain, e(a;). To derive an expression for g{uj), we ex- 
tend the procedure by Feng et al [7] by determining the 
dynamic, effective bond stiffness from a self-consistency 
requirement, as detailed below. 

The effective medium network is subjected to a macro- 
scopic infinitesimal oscillating strain e(w) = epe*"*, de- 
forming bond nm by f„me(a;). Subsequently, replacing 
this effective medium bond with one sampled from the 
distribution P{g) gives rise to an additional, nonafline 
deformation Su{u)). The original, uniform deformation 
can be restored by applying a force 



f(c.) 



(8) 



Thus, the nonafline deformation which arose from the 
bond replacement can be expressed as 



Su{uj) = 



gsAiii^) - g{uj) +g 



(9) 



where gEM (w) is the force on a bond in the effective 
medium network in response to a unit displacement. This 
allows us to express the nonafline displacement as 



Su (uj) 



fnme(^) (g (g^) - g) 

gEM (w) -g{uj)+g' 



(10) 



The self-consistency condition requires that, when av- 
eraging over all possible bond replacements, the lo- 
cal fluctuations in the deformation field must vanish, 
{Su {uj)) — 0, leading to the following equation for g{uj), 



g£;M(w) -hg-g(w) 



P{g)dg = Q. 



(11) 



We solve this equation by first determining g'j^\^{uj) as the 
displacement in response to a unit force between nodes 
n and m, f (k) = f„m (l — e*'' '"""), by solving the net- 
work's equation of motion 

u(k) = -i5-i(k).f(k), (12) 

where the dynamical matrix of the effective medium is 
given by 

g(w)r„„ (g)r„^ n^m 

(13) 

where I is the unit tensor and ® is the external prod- 
uct. As before, we set rja — 1, and the spatial Fourier 
transform of D is given by 

i?(k)^^A,e^''--'^ 

ij 

= 'Yjii.'^yij ® T^ij (l - e**'''') + iAttlxjI (14) 



Thus, the displacement of the nm bond due to a unit 
force follows 



^5]r„„.f(k)i?-Mk) (e-''"-" 



N 

k 

2g-^(c^) 

ZN 



1) 
1) 



(15) 




Y^Tr 



' i- y (1 



(1 - e*''-'") + Att 



4^1 



(1 



where Z is the maximum coordination of the underlying 
lattice, d is the dimension of the system and N is the 
total number of nodes in the network. 

For a random bond-diluted lattice, the self-consistency 
condition (Eq. ([Tl])) can be written as 



9i'^) 



^i-g (uj) 



g£;M(w) + n- g{uj) 



ii-p) — rTVi = 0' 

gEM[i^) - g{^) 

(16) 

where /i will be set to 1. By solving this equation for 
a triangular lattice configuration, we obtain the macro- 
scopic shear modulus, G*[(jj) = g(w)\/3/4. Remarkably, 
this EMT prediction for the rheology captures the main 
features of the simulation results with reasonable quan- 
titative agreement, as shown in Figs. [Sj^d. Slow conver- 
gence of the numerical integration scheme precludes a 
high-precision solution of the EMT in the critical regime. 
However, we can obtain various interesting analytical re- 
sults by considering the large or small limits of the quan- 
tity |g('^)| /w. 



High-frequency limit When uj 3> |g(tj)|, Eq. (15) can 
be written as gEM ~ 2TTUji. Using this in the self- 
consistency equation, we find the shear modulus to be 



G*(w) 



l-p 



(17) 



This high frequency limit corresponds quantitatively 
with the numerical results, as shown in the insets of 
Figs.[3}:d. 

Low-frequency limit This limit is solvable only when 
p is not much less than pc- In this case |g(w)| 3> lo and 



Eq. ( 15 ) reduces to 



9EMi^) 



2g-H^) 



2TTiuj 

g(w) 



where A is a numerical constant, 
1 1 



Er.. 



(1 



A 



5.17. 



(18) 



(19) 
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Eq. ( 20 1 can be written in the following scaling form 



Lo\Ap\- 




Lj\Ap\ 



FIG. 6. Universal scaling of the shear modulus for the numer- 
ical results (a) and the analytical calculations (b) according 
to the scaling form G* — \Ap\^Q±{Lu\Ap\~'^") over a range of 
bond probabilities. The main plots show the results for the 
elastic modulus, the insets for the loss modulus. The scaling 
exponents for the EMT calculation are / = 1 and = 2. 
An excellent collapse is observed for the numerical data with 
scaling exponents / = 1.4 and (p^ = 3.6 and Pc = 0.649 for 
a 100 X 100 network [16]. A frequency range of frequencies 
between 3.1 x lO"'' and 4.0 x 10"^ was used for the numer- 
ics, and in the EMT scaling the frequency range is 10"'' to 
{2A-Ky^ for G' and 10"^ to (2^7r)~^ for G". 



By solving the self-consistency equation for a bond- 



diluted lattice (Eq. (16l) in this limit, we find the shear 
modulus, 



6p - 4 - SiA-KUJ 



(20) 



This expression for the dynamic shear modulus captures 
the low- frequency rheology for p }^ Pc, as shown in the 
insets of Figs. [3]cd. 

The EMT indicates a critical bond probability at pc = 
2/3. Close to the critical point, |Ap| ^ 1, and in the 
limit of small frequencies {oj <C (Att)" 10~^ for G' 
and Lu ^ [Att)^ ~ 10^^ for G"), the shear modulus in 



G* = \Ap\fgi 



\Ap\' 



(21) 



where the EMT predicts f — 1 and (p^ = 2, consistent 
with the mean field predictions in ref. [551 HZ] • The scal- 
ing function Gj. {x) = Q'^. (x) -I- iQ'j. {x) is given by 



S± (x) = < cos 



„ 3\/3 . 

G± (x) = sm 



tan ^ (ax) 
tan~^ (ax) 



1 + {axY 



1/4 



± 1 



1 + (ax) 



1/4 



with a 



(22) 



branch 



g . The scaling function Q'^ has a 
above Pc and a — branch below pc- When a; — ^ and 
P > Pc G'+{x) must be constant such that G' scales as 
|Ap|-'', while G'lix) ~ x so that G" scales as u}\Ap\f-'f''-' . 
Furthermore, when a; — >■ and p < Pc GLix) ^ x^ 



such that G' scales as cij|Ap| 
such that G" 



and g'!_{x) ~ X 
scales as uj\Ap\^~'^'^ . In the critical con- 
nectivity regime, the shear modulus is finite and, thus, 

G± (x) ^ x'*'^ such that G* is Ap-independent. Conse- 
quently, G* - [iu:)^ with A = f/cj)^. 

The EMT and numerical results are collapsed accord- 
ing to this scaling form, as shown in Fig. |6] We find an 
excellent collapse for the numerical data with / = 1.4, 
determined in prior work on elastic networks [H I16j , and 
= 3.6, determined from collapsing the nonaffinity 
data (Fig. |4|. Within the EMT, A = 1/2, consistent 
with a mean field calculation in refs. [5S1 [SJ. By con- 
trast, from the collapse of the simulation data, we find 
A w 0.41. 

The scaling behavior of G* is clearly related to that of 
the nonaffinity parameter F: in both cases, the critical 
relaxation frequency is controlled by the exponent <f>^. In 
the first case, we inferred that (p^j — //A, while in the 
second case, we found that = \/5, and thus. 



/ 
A 



A 



(23) 



Solving for A and recalling that 5 = \ — A (Eq. [?]), we 
obtain 



/ 



A + / 



(24) 



Strikingly, we find that the dynamical behavior of the 
network, captured by the exponent A, can in fact be 
inferred from the rigidity exponent / and the nonaffin- 
ity exponent A of the elastic network in the absence of 
viscous interactions. Using previously determined / = 
1.4±0.1 and A = 2.2±0.4 [16], we expect A = 0.39±0.08, 
in agreement with our numerical observations. Further- 



more, using Eq. (24 1, we can also recover the mean-field 
prediction for the nonaffinity exponent A = 1, using 
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A = 1/2 and / = 1 from the EMT calculation [TU]. This 
argument should be valid for broader classes of disordered 
networks, e.g. fiber networks for which = 3.2 ±0.4 and 
Af, = 1.8 ± 0.3 in 2D, implying a dynamical scaling of 
Ab = 0.64 ± 0.13. Thus, the dynamic EMT constructed 
here, combined with the scaling arguments, provide an 
avenue for exploring dynamical behavior of a wide range 
of disordered networks. 

We thank Cliff Brangwynne for helpful discussions and 
his hospitality and Brian Tighe, Fred MacKintosh and 
Boris Shklovskii for insightful discussions. This work 
was supported by the Lewis-Sigler fellowship (CPB) and 
FOM/NWO (MS). 



[1 
[2 
[3 

[4; 

[5 
[6 

[7; 

[9 
[10 

[11 



* |cbroed er@princ eton.edu| 

" M. F. Thorpe, J. Non-Cryst. Solids 57, 355 (1983). 
J. C. Phillips, J. Non-Cryst. Solids 43, 37-77 (1981). 
H. He and M. F. Thorpe, Phys. Rev. Lett. 54, 2107-2110 
(1985). 

L. M. Schwartz, S. Feng, M. F. Thorpe and P. N. Sen, 
Phys. Rev. B 32, 4607-4617 (1985). 

S. Feng and P. N. Sen, Phys. Rev. Lett., 52, 216-219 
(1984). 

S. Arbabi and M. Sahimi, Phys. Rev. B 47, 695-702 
(1993). 

S. Feng, M. F. Thorpe and E. J. Garboczi, Phys. Rev. B 
31, 276 (1985) 

A. J. Liu and S. R. Nagel, Nature 396, 21 (1998). 
C. S. O'Hern, L. E. Silbert, A. J. Liu and S. R. Nagel, 
Phys. Rev. E 68, 011306 (2003). 

M. Wyart, H. Liang, A. Kabla and L. Mahadevan, Phys. 
Rev. Lett. 101, 215501 (2008). 

W. G. Ellenbroek, Z. Zeravcic, W. van Saarloos and M. 



[12 
[13 
[14 
[15 
[16 
[17 

[18 
[19 

[20 
[21 

[22 

[23; 

[24 

[25 
[26 
[27 
[28^ 

[29 

[30 

[31 



van Hecke, Eur. Phys. Lett. 87, 34004 (2009). 

D. A. Head, A. J. Levine and F. C. MacKintosh, Phys. 
Rev. Lett. 91, 108102 (2003). 

J. Wilhelm and E. Frey, Phys. Rev. Lett. 91, 108103 
(2003). 

C. Heussinger and E. Frey, Phys. Rev. Lett. 97, 105501 
(2006). 

M. Das, F. C. MacKintosh and A. J. Levine, Phys. Rev. 
Lett. 99, 038101 (2007). 

C. P. Broedersz, X. Mao, T. C. Lubensky and E.G. 

MacKintosh, Nat. Phys. 7, 983 (2011). 

X. Mao, O. StenuU, and T. C. Lubensky, 

arXiv:1111.1751vl (2011). 

J. G. Maxwell, Philos. Mag. 27, 294 (1864). 

M. Das, D. A. Quint, J. M. Schwarz, PLoS ONE 7(5): 

e35939 (2012). 

S. Alexander, Phys. Rep. 296, 65 (1998). 

M. Sheinman, G. P. Broedersz and F. C. MacKintosh, 

Phys. Rev. Lett. 109, 238101 (2012). 

G. Heussinger and J-L. Barrat, Phys. Rev. Lett. 102, 

218303 (2009). 

E. Lerner, G. Diiring and M. Wyart, Euro. Phys. Lett. , 
99, 58003, (2012). 

E. Lerner, G. Diiring and M. Wyart, PNAS, 109 (13), 
4798-4803 (2012). 

B. P. Tighe, Phys. Rev. Lett. 107, 158303 (2011). 

B. P. Tighe, Phys. Rev. Lett. 109, 168303 (2012). 

B. P. Tighe, arXiv:1205.2960 (2012). 

B. Andreotti, J-L. Barrat, and G. Heussinger , Phys. Rev. 

Lett. 109, 105901 (2012). 

M. Sheinman, G. P. Broedersz and F. G. MacKintosh, 
Phys Rev E, 85: 021801 (2012). 

X. Mao, O. StenuU, and T. G. Lubensky, ' arXiv: 130 1 . 08"7 
(2013). 

A. W. Lees and S. F. Edwards, J. Phys. G: Solid State 
Phys. 5, 1921 (1972). 



